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We report on a study of the classical field theory description of charged skyrmions in quantum Hall 
ferromagnets. The appropriate field theory is a non-linear a model generalized to include Coulomb 
and Zeeman interaction terms. We have tested the range of validity of the classical field theory by 
comparing with microscopic descriptions of the single skyrmion state based on the Hartree-Fock 
approximation, exact diagonalization calculations, and many-body trial wavefunctions. We find 
that the field theory description is accurate for skyrmions with moderate spin quantum numbers 
(> 10) although, as expected, it fails qualitatively for small spin quantum numbers. 



I. INTRODUCTION 

The quantum Hall cffcctEh§ occurs in a two- 
dimensional electron system (2DES) in a strong perpen- 
dicular magnetic field (B) and is associated with the ex- 
istence of incompressible ground states at certain values 
of the Landau level filling factor, i/. {v = N/N^ where 
is the Landau level degeneracy and N is the number 
of electrons in the system.) Some incompressible ground 
states are strong ferromagnets, i.e., they are total spin 
eigenstates with 5* = N/2 so that an infinitesimal Zee- 
man coupling of the magnetic field to the spin degree 
of freedom is sufficient to produce complete spin align- 
ment. Because of the small ratio of the Zeeman-splitting 
to other relevant energy scales in typical 2DES's, it is 
then useful to regard the system as a ferromagnet in a 
weak symmetry breaking magnetic field, even though the 
2DES is often in the strong field quantum limit as far as 
orbital degrees of freedom are concerned. Quantum Hall 
ferromagnets (QHF's) have a number of unusual proper- 
ties that spring from the disparity between the magnetic 
field coupling strengths for orbital and spin degrees of 
freedom. 

Because of the gap for charged excitations in incom- 
pressible states, the only low-lying excitations in these 
quantum Hall ferromagnets are those associated with 
slow variations in the unit vector field which describes 
the local orientation of the spin magnetic moment. In ad- 
dition to the spin waves modes, thece exist higher energy 
topologically non-trivial skyrmionDfl textures in the spin 
field. The importance of these topologically non-trivial 
excitations is magnified by the fact that they carry an 
electrical charge .BID As the lowest energy charge carriers 
they control the thermally activated dissipation on the 
y = 1 quantum Hall plateau and, away from fiUiiie-fac- 
tor = 1, they are present in the ground state. □ 113 In 
the latter case the ground state is no longer fully spin po- 



larized. For weak Zeeman coupling skyrmions are large 
and each charge introduced into the system flips over a 
large number of spins. Thus the spin magnetization is 
expected to show a sharp cusp at filling factor v — 1. 

The presence of these objects has recently been de- 
tected directly in NMR measurements of the electron 
spinjnagnetization which is proportional to, the Knight 
shifttJ, in thernml. transport measurementst3, and in op- 
tical absorption.EiThey may also be responsible for the 
recently observedlij enormous enhancement of the ap- 
parent specific heat because of their equilibrating effect 
on the nuclear spins. There has recently been consid- 
erable interest_in_tlie-£)hysics of these exotic electronic 
quasiparticles.B~E3ll3"L3 

Charged spin texture excitations in QHF's are the sub- 
ject of this paper. A realistic theory of these excita- 
tions requires at a minimum that Coulomb energy and 
Zeeman energy terms be added to the non-linear sigma 
(NLcr) model classical field theory. In Sec. |^ we briefly 
review the resulting classical field theory of spin textures 
in QHF's and report on numerical solutions for the low- 
est energy topologically non-trivial textures. With typ- 
ical parameters, the size of the lowest energy textures 
need not be large compared to microscopic lengths so 
that both higher order gradient corrections and quan- 
tum fluctuations not included in the model may become 
important. In Sec. Ill we present an overview of the 
various microscopic approaches which have been used to 
study skyrmion states in quantum Hall ferromagnets. In 
the Hartree-Fock approximation, the gradient expansion 
of the fleld theory is effectively summed to inflnite or- 
der. However, it becomes difficult to solve these equa- 
tions accurately for very large skyrmions. Many-particle 
variational wavefunction and exact diagonalization ap- 
proaches also incorporate, in addition, quantum fluctu- 
ation effects, but are even more limited in the size of 
quasiparticles for which accurate calculations are possi- 
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ble. Numerical comparisons of energies and spin densi- 
ties obtained using these different approximation strate- 
gies are compared in Sec. [V. One conclusion from this 
work is that the most important corrections to the min- 
imal generalized NLcr model are the higher order gra- 
dient corrections already captured in the Hartree-Fock 
approximation. The quantum fluctuations added using 
trial wavefunction or exact diagonalization are responsi- 
ble only for minor further modifications. We summarize 
our conclusions in Section |^. 



II. CLASSICAL FIELD THEORY 

A. General considerations 

Here we briefly review the classical field theory for 
charged spin texture excitations or skyrmions in quan- 
tum Hall ferromagnets.ElQ Following Sondhi et al.p we 
start from a modified version of the NLcr model clas- 
sical field theory which has been exploited to describe 
the low energy properties of two-diroepainnal Heisenberg 
ferromagnets and antiferromagnetsBjo The order pa- 
rameter in this theory is a unit vector field, m(r), thai 
describes the local orientation of the spin or pseudospinEI 
magnetic order. For quantum Hall ferromagnets the en- 
ergy functional of the minimal theory is 



E[m.] = Eo[m] + E,[m] + E,[m] 



(1) 



where -Eo[iii] is the leading order term in a gradient ex- 
pansion of the energy functional and Ez [m] is the Zeeman 
energy functional: 



Eo[m]^^ J dV(Vm)2, 



Ez[m] = 



(fr [ 1 - 7712 (r) 



(2a) 



(2b) 



Here ps — / {Idy/zneio) is the spin stiffness (assuming 
zero layer thickness for the 2DEG) which can be cal- 
culated analytically for this case and represents a loss of 
Coulomb exchange energy, t = {g*fiBB)/2 represents the 
Zeeman coupling strength, e is the dielectric constant of 
the host semiconductor, and £o is the magnetic length. 

When only these first two terms are present, the spa- 
tial extent of structure in skyrmion extrema of the energy 
functional shrinks to a point and the gradient expansion 
fails. To obtain a consistent theory of skyrmion excita- 
tions when the Zeeman coupling is non-zero it is neces- 
sary to go beyond leading order in the gradient expan- 
sion. For QHF's it turns outcl that, because skyrmions 
are charged, the next to leading term in the gradient 
expansion is the non-local Coulomb interaction energy 
functional 



, P(r)p(i-') 



(2c) 



where the charge density is given b 



p(r) = eab m(r) • [9am(r) A 9f,m(r) 



(3) 



The charge density is the product of the Landau level 
filling factor of the ground state v and the 0(3) topo- 
logical density of the spin texture. The non-locality of 
this functional complicates the calculations we present 
below. The total skyrmion charge is therefore v times an 
integer- valued topological invariant. If we compactify the 
2D plane to a sphere 5^, then m(r) defines a mapping 
of S'^ onto the order parameter sphere S'^ and hence can 
be classified by the homotopy group 7r2(S'^). The topo- 
logical charge simply counts the wrapping^joLlhe order 
parameter sphere by the coordinate sphereHj'EJ 

Extrema of the energy functional of this model, m, 
satisfy a non-linear differential equation which can be 
obtained by minimizing Eq. (|l|) with respect to m, using 
a Lagrange multiplier to enforce the constraint m(r) • 
m(r) — 1. After eliminating the Lagrange multiplier we 
find that 



-V^ + m • V^m)m^ - 



27r£( 



m^m,, - 



47r 



eafcRnr)}(mAa,m)^ = 0, (4) 



where V{y) is Hartree potential 



V(v) = ^ / 
eo J 



dr 



, P(r') 



(5) 



and p is the skyrmion charge density corresponding to the 
minimum energy solution, m(r). The solutions of Eq. (0) 
can be classified by the skyrmion charge Q = J dr p{r). 
In the absence of Zeeman and Coulomb energies the en- 
ergy functional is scale invariant and Eq. (Ht-b?^ a family 
of known analytic scale invariant solutions!tll"E3 The so- 
lutions may be represented by analytic complex valued 
polynomials: 



m(z) 



4wi 



4:W2 



4' 



(6) 



where w — Wi + iw2 is a polynomial in the complex 
variable 2 = (x + iy)/\ with arbitrary scale A. The de- 
gree of the polynomial w determines the skyrmion charge. 
These analytic solutions are not valid if the Zeeman and 
Coulomb terms are included. In the following subsection 
we briefiy explain the numerical methods we use to solve 
the generalized differential equation for the case of unit 
charge spin texture. 



B. Single skyrmion 

In this subsection we concentrate our attention on 
antiskyrmions (skyrmions) with unit topological charge, 
Q = 1(— 1). We find a numerical solution to Eq. (^ for 
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an infinite system with the boundary condition that the 
spin is down at the skyrmion center and up at infinity. 
It is convenient to choose the origin of coordinates at the 
center of the skyrmion to take advantage of the circular 
symmetry of the charge distribution of a single skyrmion. 
We parameterize the order parameter m(r) by spherical 
angles so that 

m(r) = (sin6'(r) cos (/)(r), sin 0(r) sin (/>(r), cos 0(r)) , (7) 

where d{r) and 0(r) are new field variables. We take ad- 
vantage of the circular symmetry for the single-skyrmion 
problem by choosing 0{r, ip) = Oir) and 4>{r, ip) = zLip+ipQ 
where ip is the azimuthal angle and ipo is a constant. 
The energy of the skyrmion is independent of ipo be- 
cause of the invariance of the energy functional under 
global rotation of all spins in the plane about the z- 
axis. Positive and negative signs for p yield antiskyrmion 
and skyrmion solutions, respectively. One may show 
explicitly that the above ansatz and boundary condi- 
tions yield unit winding number, Q = ±1, since the 
charge density in the circularly symmetric case is given by 
p{r) = ^ (1 / Anr) {d / dr cos 9 (r)). It follows that skyrmion 
and antiskyrmion solutions have the same energjj, allow- 
ing us to consider skyrmion spin textures only.Eil In this 
case, Eq. reduces to a non-linear integro-diffcrential 
equation for 0{r): 

loTT^ r J dr dr 

Here 

U{ry)^- -K[} M (9) 

where K{x) is the complete Elliptic integral of the first 
kind continued appropriately to negative values of x. 
We have solved Eq. (||) for 9{r) using an iterative ap- 
proach. At each step in the calculation the integral over 
r' was evaluated at each r using an approximation for 
9{r). When this integral is fixed, we are left with a 
two-point boundary value problem [9(r = 0) = tt and 
9{r oo) = 0] which can be solved by standard meth- 
ods. The resulting value of 9{r) to evaluate the integral 
over r' for the next iteration. We found that this it- 
erative procedure converged rapidly when the iteration 
was started from one of the analytic solutions obtained 
neglecting Coulomb and Zeeman energies. 

Figure |l] illustrates some of the results obtained from 
these numerical calculations. Wc plot the size of the 
skyrmion (A), defined by 9[\) ~ n/2, and the the number 
of reversed spins (K), defined by 

K^'Jdr[l-m.ir)]-l (10) 



as a function of t. We find that A ~ t~^/^ and K 

The logarithmic corrections to these power laws predicted 

on the basis of a variational solution of the field theorycl 

are not yet apparent at the smallest values of t we have 

considered. 
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FIG. 1. Skyrmion size (circles) and spin (black dots) ob- 
tained from the classical field theory as a function of the Zee- 
man coupling constant t. The skyrmion size A is defined as the 
radius at which the spin lies in the xy plane, i.e., cos[6'(A)] = 0. 
These results are consistent with the power law oc t~^^^ be- 
havior expect at small t. The number of reversed spins in a 
Skyrmion is proportional to A^ cx: t~^^^. Its deviation from 
these power law at large t is an indication of the necessity for 
including microscopic physics not captured by the minimal 
field theory. 

Figure || presents our numerical results for the depen- 
dence of skyrmion energy on Zeeman coupling strength 
in the minimal field theoretical description. The Blavin- 
Polyakov solutions of Eq. (^ have the minimum possible 
gradient energy, which has the value When the 

Zeeman and the Coulomb terms are added the skyrmion 
texture is reshaped and the optinuil solutions are no 
longer analytic in z. In particulam the component of 
the magnetization perpendicular to the Zeeman field has 
an exponential rather than a logarithmic fall of at large 
distances. These non-analytic solutions must have a 
gradient energy larger than 4,71 ps, with the deviation 
from this value decreasing as the Zeeman energy de- 
creases, or, equivalently, as the size of the skyrmion in- 
creases. We thus expect that for the gradient energy 
Eo{X) — 'inps + /(A) where /(A) is monotonically de- 
creasing function of the skyrmion's size (see Fig. ^) . The 
non-analytical solutions become asymptotically close to 
the Blavin-Polyakov solutions at small Zeeman splitting 
(large A). Dimensional analysis shows that, up to possi- 
ble logarithmic factors, the Coulomb and the Zeeman en- 
ergies must vary at large A as C/A and DX^, respectively, 
where C oc e^/e and D oc t/ig are constants. Assuming 
that /(A) vanishes faster than 1/A, optimizing the total 
energy [47rps + /(A) + C/A + DA^] with respect to A leads 
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to the prediction = Ec/2 oc {e^ / ee.of'^t^/^ as A -> oo. 
These power laws were originally suggestedQ by Sondhi 
et al. Our numerical results confirm these predictions, 
which hold with surprising accuracy out to fairly large 
values of t: in fact E^/Ec is within 1% of its asymp- 
totic value in the whole range of values of A we have 
studied. The fact that /(A) vanishes faster than A~^ is 
expected on the basis of perturbative treatments of the 
Coulomb and Zeeman terms and confirmed by the numer- 
ical results shown in Figure |. The total energy of the 
Skyrmion is closely given by E/pg = 47r + yl[g| ln((7)|]-'^/'^, 
where g — 2t/(e^/e£o) and, from Figure ||, we estimate 
the dimensionless constant A « 30. (Unlike the results 
for A and K , the logarithmic factor in the energy is ap- 
parent in our numerical results at small Zeeman coupling 
strengths.) This numerical estimate|.aL|j4 compares rea- 
sonably well with previous estimatesQo in which the re- 
sult A « 25 is obtained. The consistency of these numer- 
ical results affords confidence in their accuracy. 



-* tog[E/p^ 
log[£/p, 
log[£'/p^-47c] 



1.0 



0.0 



-1.0 




-5.0 



-4.0 



-3.0 -2.0 

iogk|inte)|] 



-1.0 



0.0 



FIG. 2. Skyrmion energies: Ec (crosses), Ez (stars), ex- 
cess gradient energy Eq — Att (circles), and total excess energy 
(black dots) as a function of log[(7ln((;)] [g is the Zeeman spin 
splitting in units of e^/(e^o)]. The power law for the depen- 
dence of the excess skyrmion energy, which we extract from 
a linear fit to these plots, is close to 1/3. 



the direction of the field. All the eigenstates of the sys- 
tem can be labeled by the quantum numbers S and S^. 
As mentioned in the Introduction, the particular case 
N — corresponds to ^ = 1 and, although there is no 
rigorous proof, numerical evidence clearly indicates that 
the ground state has the maximum possible total spin 
S = N/2 and, in the absence of a Zeeman field_a degen- 
eracy 25* -I- 1. As pointed out by Jain and WuE3, this is 
one of the few cases where the first Hund's rule applies 
in the fractional Hall regime. When the Zeeman field is 
in the z direction, the ground state has Sz = N/2 and 
the electronic system is completely spin-polarized in the 
direction of the field at zero temperature. 
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FIG. 3. Energy spectrum obtained from an exact diago- 
nalization of the Hamiltonian in the spherical geometry with 
A'' = 8 and A^,^ = 9 for (a) hard-core and (b) long-range 
(Coulomb) interactions (no Zeeman coupling included). Dif- 
ferent symbols have been used for different values of \Sz\. 
For the case of the hard-core model a degenerate band cor- 
responding to the Landau level degeneracy occurs once for 
each member of each S — N/2 — K multiplet present at this 
system size. For the case of long-range interactions the de- 
generacy is removed and the energy decreases as K{S) in- 
creases(decreases) . Rotational symmetry is always present in 
the spherical geometry and the Landau level degeneracy of 
the charged quasiparticles is preserved for each value of K. 



III. MICROSCOPIC DESCRIPTIONS OF SINGLE 
SKYRMION STATES 

A. Exact diagonalizations 

The microscopic physics of skyrmions can be ap- 
proached directly by numerically diagonalizing the finite- 
size Hamiltonian matrix for N two-dimensiona|l S L pij i-1/2 
electrons in A^^ lowest Landau level orbitals.crEjO The 
Hamiltonian commutes with the total spin operator, S"^, 
and, when the Zeeman term is neglected, also with its 
projection, S^, on any direction ^. When the Zeeman 
coupling is present the direction ^ must correspond with 



Since skyrmions are charged quasiparticles, they ap- 
pear in exact diagonalization calculations performed with 
A^ — Nij, ± 1. (Particle- hole symmetry permits exami- 
nation of one sign of charge only.) In order to-obtain 
the energy spectrum in the spherical geometryl3l23'EZI we 
need all interaction iftatrix elements in the lowest shell of 
monopole harmonics" which comprise the lowest Landau 
level basis. These are given by 

{mim2\V{n - f1')l"i3m4) - ^(-l)™+"H+™2 

y ( ' ' (11) 

^-^ \ ^ra\ — m 7714 / \ — m2 m 

/=|m| ^ ^ ^ ^ 
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where m = — mi = TO2 — and the rrii are the 
azimuthal angular momentum quantum numbers which 
label the single particle orbitals, V{(l — (l') is the inter- 
action potential where fl denotes the electronic spherical 
coordinate, and Ns = {N^ — l)/2. The coefficients in 
this expansion are products of Wigner's 3-j symbols and 
the Vj^^ parameters are the interaction-dependent pseu- 
dopotentials introduced by Haldane. Since the Hamil- 
tonian also commutes with the total orbital angular mo- 
mentum L the diagonalization is carried out in subspaces 
of fixed Sz and M = L^- It is instructive to consider 
first the hard-core model for the electron-electron inter- 
actions. In this case the pseudopotentials V/^^ are given 

by 



(2iV5 + 1)2(2/ + 1) 



Ns I Ns 
-Ns Ns 



(12) 



For this model, as illustrated in Fig. ^(a) and discussed 
previously by Jain and Wulj and MacDonald, Fertig, and 
Breycl, multiplets of zero energy eigenstates occurs with 
all possible values of the total spin quantum number S. 
In a quantum description, the number of reversed spina 
in a skyrmion spin texture is quantized and conjugateO 
to its global orientation. The zero-eMrgy spin-multiplets 
with S = N/2 — K have been shownEI to be the quantum 
states which correspond to a classical skyrmion texture 
with K reversed spins. The results shown in Fig. ^(b) are 
for the case of Coulomb interactions between the spins, 
for which the pseudopotentials are 



V, 



Ns 



{2Ns 



Ns 

-Ns 



Ns 
Ns 



(13) 



In this case skyrmion states with more reversed spins 
have lower energy, at least when the Zeeman energy is 
neglected. When the Zeeman energy is included, the en- 
ergy will in general be minimized at an intermediate value 
of K. This is the result first obtained by RezayiE^, which 
provided the first evidence for the existence of charged 
excitations with large spin quantum numbers in quantum 
Hall ferromagnets. The appearance of degenerate states 
with different total angular momenta reflects the exten- 
sive Landau level degeneracy of single-skyrmion states 
with all values of K. 

We have also performed the same calculation {N = 8, 
N^ = 9) in the disk geometry. Figure ^(a) shows the 
low-energy spectrum for the hard-core model where one 
can see the same degenerate band at zero energy as in 
the spherical geometry case. For the zero-energy states 
of the hard-core model the Landau level degeneracy is 
preserved. For the Coulomb interaction, on the other 
hand, Landau level degeneracy is broken by strong edge 
effects [see Fig. §(b)] and the non-local nature of the 
interaction. For a given value of M , the minimum in- 
teraction energy state is always found at the maximum 
allowed value of K. 
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FIG. 4. Energy spectrum obtained from an exact diago- 
nalization of the Hamiltonian for a disk geometry with N = 8 
and = 9, for (a) hard-core and (b) long-range (Coulomb) 
interactions (no Zeeman coupling included). Different sym- 
bols have been used for different values of 1. As in Fig. ^ 
for the case of the hard-core model, a degenerate band at zero 
energy is observable where all possible values of K are present. 
Unlike in the spherical geometry for the case of long-range in- 
teractions the translational symmetry is broken due to strong 
edge effects. 



B. Hartree-Fock 

A simpler and perhaps physically more transparent mi- 
croscopic approach for computing skyrmion properties is 
the unrestricted Hartree-Fock approximation approach. 
This procedure finds the lowest energy Slater determi- 
nant wavefunction which supports a spin texture of the 
form found in skyrmion solutions of the NLcr model. Such 
wavefunctions gcnerically have the formEl (for quasihole 
skyrmions) 



I* 



HF) — 



t 



it)|0), 



(14) 



where cj,j^ creates a lowest Landau level electron in an- 
gular momentum state m with spin a, and the Um, Wm's 
are variational parameters obeying the constraint |umP + 
= 1. The precise way in which the energy of 
the state I^'hf) may be minimized with respect to the 
Um, Vm's has bccu discussed elsewhereO. 

In practical calculations, the Hartree-Fock approach 
bridges the gap between skyrmion sizes which can be 
reached with exact diagonalization approaches and sys- 
tems sizes where the field theoretic approach becomes 
accurate. Typically one is limited to system sizes with 
N^ « 10 in finite size diagonalization studies, since the 
matrix dimensions grow rapidly with system size. On 
the other hand, Hartree-Fock calculations are possible for 
spin textures involving up to several thousand electrons. 
As in exact diagonalization studies one works with the 
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exact many-body Hamiltonian in Hartree-Fock studies 
so that the quantitative errors introduced in the gradi- 
ent expansion of the energy functional used to gengrate 
the NLfT model for the quantum Hall ferromagnetHlJ are 
not present. The Hartrcc-Fock approach also captures 
some of the quantum mechanics of the problem, though 
clearly not all. In particular, I'^'hf) has M — Sz, the 
difference of the orbital and the z-component of the spin 
angular momenta, as a good quantum number. How- 
ever, M and Sz are not separately quantized as they 
should be, nor is |^h_f) an eigenfunction of 5^. These 
failures could easily have some quantitative importance 
for K ^ 3, the quantum number of the lowest energy 
skyrmion in typical experiments. From a field theoretic 
point of view, one may regard the Hartree-Fock approx- 
imation as a mean field theory which retains the higher 
order gradient terms absent in the simpler field theory 
studied in the previous section. The resulting effective 
action would contain most information, but, in principle, 
one should include quantum fluctuations around this sta- 
tionary phase approximation in order to properly account 
for the quantization of orbital and spin angular momen- 
tum. As shall be seen below, for the purposes of comput- 
ing energies and sizes of the skyrmions, the Hartree-Fock 
approach appears to be remarkably accurate, suggesting 
that such quantum fluctuations are quantitatively unim- 
portant. However, for probes sensitive to the quantum 
numbers of the skyrmions - tunneling spectroscopy for 
example - wavefunctions with the correct quantum nupi- 
bers are needed to obtain qualitatively correct resultsc2l. 



C. Variational wavefunctions 

An alternate quantum description of single-skyrmion 
states can be obtained by using variational many-body 
wavefunctions. These can be refined by comparing with 
numerically exact solutions for small K and can provide 
a convenient, if somewhat indirect, route to calculate 
properties of physical interest — » oo limit. The vari- 
ational wavefunctions we employ are motivated by the 
disk-geometry wavefunctions introduced by MacDonald, 
Fertig and Breya which, in the thermodynamic limit, 
become exact for hard-core-interaction single-skyrmion 
states. The second quantized form of these wavefunctions 
clearly exhibits the microscopic nature of the skyrmion: 



\K) = 



1 



E 



I 



\/C(K) , Jnii ■ ■ ■ niK 

cLk-uC™a-T • • -cL^-uC^iTCotI*)- (15) 

In the above expression \^) denotes the ferromagnetic 
11 = 1 state and C'{K) are the normalization constants 



C{K) = 



mK>->nii = l 



mi ■ ■ ■ niK 



(16) 



This wavefunction correspond to a skyrmion with quan- 
tum number K located at the origin of the plane with 
Mmax — Mfioie — K [Mhoie corrcsponds to the total an- 
gular momentum with a bare quasihole located at the 
origin [see Fig. ^(a)]). All other single-skyrmion states 
with different orbital angular momenta M within a par- 
ticular zero-energy K band can be obtained by lowering 
the center-of-mass angular momentum n 
times: 



Kl 



(17) 



In a large system where edges may be neglected, all such 
states are degenerate. Here we focus only on the states 
of highest orbital angular momentum, which are centered 
at the origin. 

It is interesting to note that these wavefunctions are 
identical to hard-core model Hartree-Fock wavefunctions 
[Eq. (^)] projected onto a state of definite S^. (For 
the hard-core model, the coefhcients Um,Vm m Eq. jl^ ) 
may be determined analyticallyij.) Thus, they presum- 
ably improve on the Hartree-Fock wavefunction in that 
they include the quantum fluctuations necessary to ob- 
tain the quantization of spin. It should also be noted 
that Eq. ( [l5|) is not precisely an eigenstate of 5^ , as may 
be verified by acting on it with the total spin raising op- 
erator S^. This operator should but fails to annihilate 
the state for any finite size system. However, it is eas- 
ily seen that S'+j^'^) cx C(if)~^/^, which vanishes in 
the thermodynamic limit. Thus, the relative weight of 
this wavefunction outside its appropiidate spin multiplet 
vanishes with increasing system size.c3 



-0.10 



-0.15 



-0.25 
0.00 



• • • 




•k=o* 



0.02 0.04 0.06 0.08 

1W. 



K=3 • 

0.10 



FIG. 5. Energy gain of the skyrmion as a function of 
for different values of K. These results were obtained using 
unmodified hard-core model wavefunctions in the disk geom- 
etry (see text) . 

We use these variational wavefunctions to estimate the 
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energies of single skyrmion states for the case of Coulomb 
electron-electron interactions. Because of numerical dif- 
ficulties created by the long range of this interaction it 
is convenient to introduce the neutral excitation energy 
defined as follows: 



e(i^) = e+(0)+e_(X) 



(18) 



where e-i-(O) is the energy to make a K = antiskyrmion 
(bare spin- reversed electron), 

e+(0) = E{K = 0,N = N^ + 1)~ E{N = N^), (19) 

and e-{K) is the energy to make an skyrmion with K 
reversed spins, 

e^iK)^EiK,N^N^-l)-EiN = N^). (20) 

The quantity of physical relevance is the difference e{K) — 
e(0), the interaction energy gain resulting from a texture 
with K reversed spins. The calculation of the neutral 
excitation energies has been done for up to = 32 and 
the results for this energy gain are shown in Fig. ^ as a 
fimction of The edge effects, which bear primary 

responsibility for the finite-size dependence seen in Fig. 
^, are minimized by increasing the angular momentum 
of the minority-spin electron in the K = 0, N = N^p + 1 
state in step with K. This allows us to extract reliable 
results from calculations at fairly small values of .cil In 
order to make an accurate extrapolation of these results 
to the thermodynamic limit we use a Wynn's e algorithm 
to generate Fade approximantscj The results of this ex- 
trapolation are shown in Fig. |^, where the same method 
has been used to extrapolate the exact resul ts ob tained 
in spherical geometry calculations (see Sec. [II A ) with 
Nff, as large as 11. 
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FIG. 6. Energy gain of the skyrmion as a function of K. 
The plot shows energies obtained using the hard-core wave- 
functions, the optimized variational wavefunctions, and en- 
ergies obtained from exact diagonalization in the spherical 
geometry. We were able to reliably extrapolate the exact di- 
agonalization results only for K up to 3. 



The substantial differences between exact diagonaliza- 
tion and variational wavefunction results is indicative 
with deficiencies in these wavefunctions. The source of 
the difficulty is that these wavefunctions do not repre- 
sent a finite size object as iV^ oo. (This can be eas- 
ily seen from the fact that the normalization constants 
C{K) diverge in the thermodynamic limit.) From coef- 
ficients in the expression ([l^) one can see that the tail 
in the spin density distribution falls off like l/r^, which 
coincides with the long-distance behavior of the classi- 
cal soliton solution in the conventional NLa model for 
fcrromagnets. However, as stated in Sec. |l[ the long- 
distance behavior departs from this form when Coulomb 
and Zeeman terms are present. In the classical solution 
the 1/r^ behavior is replaced by an exponential decay at 
large distances. This observation motivates the following 
refinement in our variational wavefunctions: 



m(«)> 



E 

:>->m 



g-Q(r?JiH hmjf ) 



\/A(K) , ^TOi • • • niK 



-uCm^T ■ • •cLi_i|C™itCoT|*>, (21) 



where the A(K) are the new normalization constants 



mjc>..>mi=l 



-2a(miH hrriA') 



mi • ■ • rriK 



The parameter a offers the possibility of optimizing the 
shape of the skyrmion for each value of K. 

It is, perhaps, surprising that the introduction of a 
Coulomb repulsion leads to a more rapid fall off in the re- 
versed spin density at large distances from the skyrmion 
center. One way of understanding this difference comes 
from examining the wavefunctions in Eqs. ( p^ and ([2l|). 
For K = 1 these are sums of products of a quasihole in the 
m = orbital and a particle-hole pair with the minority- 
spin electron precisely one unit of angular momentum 
lower than majority-spin hole. Since such particle-hole 
pairs are precisely what are needed to construct spin- 
waves at v — lc3, it is convenient to think of the state 
as a spin wave with orbital angular momentum m = — 1 
moving in the presence of a spin-polarized quasihole. For 
the case of hard-core interactions, it is easy to see that 
the spin-wave disturbance and the quasihole do not inter- 
act. The spin wave is therefore unbound and the linear 
scale over which the spin disturbance is present is lim- 
ited only by the system size. For Coulomb interactions, 
the long-range nature of the interaction leads to quite 
a different result. Because the spin-minority electron is 
always closer to the origin than the spin-majority hole, 
there is an effective attractive interaction between the 
spin wave and the quasihole at the origin. The skyrmion 
thus represents a state in which the spin-wave is bound 
to the quasihole so that the spin disturbance is local- 
ized near the origin. Thus one should expect the weight 
for large momentum particle-hole pairs to fall off more 
quickly than for the case of the hard-core interaction. 
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Introduction of the parameter a provides the variational 
freedom necessary for this to occur. 

A deficiency of the states in Eq. (|2^) is that they are 
not eigenfunctions of S'^ (although they are, of course, 
eigenstates oi Sz). They can be modified to have correct 
quantum numbers at a considerable cost in complexity 
as discussed in Ref. interested readers are referred 
to that work for details. The results we obtain using 
the present wavefunctions extrapolating to the thermo- 
dynamic limit and the values of a that optimize the en- 
ergies are shown in Fig. ^; the marked improvement over 
the hard-core interaction wavefunctions (^_5|) is apparent. 
As the value of K is increased, the exponential factor be- 
comes irrelevant and one recovers the long-distance be- 
havior expected for classical solitons. This limit has been 
studied in Ref. |l6| and Ref. 15 and tJae result obtained for 
the energy of the classical skyrmionO is reproduced by the 
hard-core wavefunctions. 



IV. COMPARISON 

This section is devoted to comparison of results ob- 
tained from the classical field theory, Hartree-Fock ap- 
proximation, and variational wavefunction for Coulomb 
interactions between electrons. In Fig. ^ we plot results 
obtained for A(ii') = e{K) — e{K + 1) calculated with all 
three approaches. 
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FIG. 7. Skyrmion energy differences A{K) obtained from 
exact diagonalization (circles), variational wavefunctions 
(black dots), the Hartree-Fock approximation (dashed line), 
and the field theory (solid line). One can see the relative 
importance of including higher-order gradient contributions 
(Hartree-Fock) and quantum fluctuations (variational wave- 
functions and exact diagonalizations) . A(_ftr) is the Zeeman 
spin-splitting value at which the number of reversed spins in 
the lowest energy skyrmion changes from K to K — 1. 



fine A{K) = de/dK] A(ivr) is the value of the Zeeman 
spin splitting (2t in the notation of this paper) energy 
at which the single skyrmion states with K and K ~~ 1 
reversed spins are equal in energy. For GaAs systems 
t^{K) I ((? I tl) K, 0.006^B[Tesla], depending somewhat 
on the details of the semiconductor quantum well so that 
K is typically ^ 3 for experiments performed at fields 
^ lOTesla. It is clear from this plot that the corrections 
to the minimal field theory captured by the Hartree- 
Fock approximation are sufficient to yield accurate re- 
sults. The use of leading gradient terms in the minimal 
field theory is seriously in error for the small K states of 
practical interest. The additional quantum fluctuation 
effects included in the variational wavefunction and ex- 
act diagonalization calculations are never of overriding 
importance. 

Figure ^ compares the spin-textures by plotting vriz (j) 
for classical field theory, Hartree-Fock, and variational 
wavefunction approaches at small integer values of K. 
For the classical field theory and Hartree Fock approxi- 
mations we plot the solution of Eq. (||), which gives inte- 
ger values for the number of reversed spins K as defined 
in Eq. ([L0[). In the quantum calculations 



EE 

m s— ±1 



(22) 



where n^^ = {'^\{a)\c\^„Cry^a\^%i{oL)) are single particle 
occupation numbers, and V'm(r) are the angular momen- 
tum states of the lowest Landau level in the symmetric 
gauge. Note that for the if = skyrmion state, the vari- 
ational wavefunction and the Hartree-Fock calculations 
are identical. In fact this state is a simple spin-polarized 
quasihole state which is known to be given exactly by 
a single Slater determinant. Apparently the fact that 
the Hartree-Fock approximation becomes exact both for 
= and K ^ oo leads to good accuracy at all values of 
K. As expected, the field theory results approach those 
obtained with the variational wavefunctions and Hartree- 
Fock approximations for large enough skyrmions while 
for small skyrmions there are significant differences. The 
discrepancy is largest for K — ^ where the classical field 
theory calculation indicates that the spin is reversed at 
the origin. In fact this property was one of the boundary 
conditions used to solve the Euler differential equation of 
the classical field theory. On the other hand, the quan- 
tum calculations show that the z-component of the spin is 
zero at the origin; in fact the total spin-density is zero at 
this point because the total charge density is zero. These 
changes in the local properties are outside the scope of 
the generalized NLct model energy functional. These er- 
rors are important for the smallest skyrmions but be- 
comes less significant with increasing skyrmion size. 



[For the classical field theory and the Hartree-Fock ap- 
proximation calculations, K is not quantized and we de- 
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rli, r/Z, 



FIG. 8. Radial distribution of in three approximations: 
classical field theory (solid lines), variational wavefunctions 
(dotted lines), and Hartree-Fock (dashed lines). The four 
smallest spin textures have been considered: K — 0(a), 
K = 1(b), K = 2(c), and K = 3(d). For K = the re- 
sults of the variational wavefunctions and the Hartree-Fock 
approximation lead to the same solution and, as expected, 
both schemes converge to the classical field theory for increas- 
ing skyrmion size. 

On the other hand, the quantum exact diagonalization, 
variational wavefunetion and even Hartree-Fock calcula- 
tions cannot be carried out for very large skyrmions. The 
limit for the many-body calculations is K ^ 4 vifhich is 
adequate for comparing with existing experiments but 
not for addressing asymptotic values of various quantities 
for the very large skyrmions which would be obtained if 
the Zeeman coupling could be adjusted experimentally to 
small values. Finite-size effects in these quantum calcula- 
tions are exacerbated by the long almost algebraic tails of 
large K skyrmion states. For example, Hartree-Fock cal- 
culations with > 1000 are necessaryS to obtain a ac- 
curate values of A(ivr) near K = 50. Similarly it has not 
proven possible to directly demonstrate the asymptotic 
behavior , E — > Airps , as i — > using Hartree-Fock calcu- 
lations. (However, the asymptotic value may be demon- 
strated by subtracting the Hartrce contribution to the 
Hartree-Fock energy and noting that this should vanish 
as the skyrmion size diverges.) For the large skyrmion 
states, which hopefully will be studied in future experi- 
ments, the classical field theoretic approach is accurate 
and convenient. 

The discrepancy between the results of the classical 
field theory and the microscopic calculations could be 
reduced by including quantum fluctuations around the 
classical solution of the minimal field theory. A uniform 
ferromagnet does not have quantum fluctuations around 
its fully spin-aligned classical ground state and this is 
also the exact quantum ground state. However, near a 
skyrmion center, the spin direction changes with position 
and quantum fluctuations will occur. As the size of the 
skyrmion is reduced, these fluctuations will become more 



important. It is clear that the sign of the effect will be 
to reduce the spin polarization near the origin and bring 
the field theory results into closer agreement with the mi- 
croscopic results. SucK-calculations have been done for 
the simple NLcr model.ca However jComputation of these 
fluctuations is numerically difficultCj, even at the Gaus- 
sian level, in the present model due to the non-locality of 
the Coulomb interaction. 



V. CONCLUSIONS 

We have presented a study of single skyrmion states of 
quantum Hall ferromagnets using and comparing a vari- 
ety of techniques. We have obtained the classical saddle 
point solution of the appropriate non-linear sigma model 
modified to include the long-range Coulomb interactions 
and Zeeman coupling. This minimal classical field theory 
model includes the leading order gradients necessary to 
obtain stable skyrmions in the presence of Zeeman cou- 
pling. We have developed and evaluated variational wave 
functions describing states with the appropriate quantum 
numbers. Finally, we have performed microscopic numer- 
ical exact diagonalization calculations on finite size sys- 
tems. These results are compared with results obtained 
using a Hartree-Fock approximation and reported in ear- 
lier work. 

Our classical field theory results confirm that the 
skyrmion size scales as the inverse cube root of the 
Zeeman energy as expected from the competition be- 
tween the Zeeman energy and the Coulomb energy and 
predicted earlierQ by Sondhi et at. We find that, for 
skyrmion states with small numbers of reversed spins, 
the Hartree-Fock results are in excellent agreement with 
many-body exact diagonalization calculations and with 
calculations based on variational wavefunctions moti- 
vated by exact zero-energy eigenstates of hard-core model 
systems. Comparison of these calculations suggests that 
the minimal classical field theory model is not quantita- 
tively adequate for the size of skyrmion which plays an 
important role in typical experiments in GaAs quantum 
well systems. On the other hand quantum Hartree-Fock 
calculations are shown to be in excellent agreement with 
many-body exact diagonalization and variational wave- 
function approaches. Large skyrmion systems may be- 
come experimentally available if the Zeeman coupling 
strength can be tuned to smaller values, for example by 
application of uniaxial stress to semiconductor host of 
the quantum well, and disorder effects can be controlled. 
In this case quantum calculations become impractical 
and the classical field theory approach, possibly supple- 
mented by field theory quantum fiuctuation calculations, 
may be necessary to describe experiments. 
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